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The formalism developed previously for the calculation of the far-wing line shape function and the 
corresponding absorption coefficient using a single-component anisotropic interaction term and the 
binary collision and quasistatic approximations is generalized to multicomponent anisotropic 
potential functions. Explicit expressions are presented for several common cases, including the 
long-range dipole-dipole plus dipole-quadrupole interaction and a linear molecule interacting with 
a perturber atom. After determining the multicomponent functional representation for the interaction 
between C0 2 and Ar from previous published data, we calculate the theoretical line shape function 
and the corresponding absorption due to the v 3 band of C0 2 in the frequency region 2400-2580 
cm -1 and compare our results with previous calculations carried out using a single-component 
anisotropic interaction, and with the results obtained assuming Lorentzian line shapes. The principal 
uncertainties in the present results,* possible refinements of the theoretical formalism, and the 
applicability to other systems are discussed briefly. 
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I. INTRODUCTION 

During the past decade, substantial progress has been 
made on the longstanding problem of the accurate calcula- 
tion of the absorption of radiation by the far wings of al- 
lowed spectral lines. Since the publication of a quasistatic 
theory for the line shape calculated within the binary colli- 
sions approximation by Rosenkranz, 12 several authors have 
extended and generalized this formalism. Boulet et al? pub- 
lished a line-by-line generalization of the quasistatic theory, 
and Ma and Tipping, in a series papers, 4-12 have improved or 
eliminated many of the restrictions and assumptions that lim- 
ited the accuracy and applicability of the original Rosen- 
kranz theory. Although the theoretical results for the specific 
cases considered (H 2 0-H 2 0, H 2 0-N 2 , H 2 0-C0 2 , and 
C0 2 -Ar) are in substantial agreement with both laboratory 
data and atmospheric measurements, some discrepancies re- 
main. 

All the calculations discussed above have been carried 
out using only a single-component anisotropic interaction 
that is responsible for the broadening. Because of the sensi- 
tivity of the results to changes in the anisotropic interaction, 
both the magnitude and the temperature dependence of the 
line shapes and the absorption, and also because of the inap- 
plicability to represent well the anisotropic potential by a 
single-component anisotropic model, we present in the 
present paper a generalization of the quasistatic formalism to 
the multicomponent anisotropic potential. The new formal- 
ism is presented in Sec. II. In Sec. Ill, we discuss Several 
cases, including the long-range dipole-dipole plus dipole- 
quadrupole interaction and a linear molecule interacting with 
a perturber atom. After obtaining a representation for the 
interaction between a C0 2 molecule and an Ar atom, we 
apply the theory in Sec. IV to the high-frequency wing of the 
n 3 band of C0 2 ; we compare the present results for the ab- 


sorption coefficient in the 2400-2580 cm -1 spectral regions 
to those obtained from a similar theory using only a single- 
component anisotropic model and to those obtained by as- 
suming Lorentzian line shapes. The conclusions concerning 
the present work, the principal theoretical uncertainties and 
possible refinements of the formalism, and the applicability 
to other systems are discussed briefly in Sec. V. 

II. GENERAL FORMALISM 

A. The absorption coefficient and spectral density 

For a low-density gas sample, we divide it into absorber 
molecules and the remaining bath molecules and we focus on 
one absorber molecule only. Then, with n a absorber mol- 
ecules per unit volume, the absorption coefficient per unit 
volume can be expressed as 

4i r 2 

a ( <t) ) = Ti. ~ n a to tanh(ha>/2kT)[F{ a)) + F (- «)]. (1) 
one 

In this expression, F{u>) the spectral density is given by 

F(u)=- Re Tr f% , '“ , </4" ) (0 (2) 
tr Jo 

where the angular brackets denote the ensemble average over 
all (one absorber molecule plus bath) variables. The dipole 
operator of an absorber molecule in the Heisenberg represen- 
tation /tio W) (r) is determined by 

M W(t) = e iH,/A M i H) (0)e- iH,/A = e iL, M l a H) (0). (3) 

In the above expression, H is the total Hamiltonian which 
consists of the unperturbed Hamiltonian of the absorber mol- 
ecule and the bath, H 0 = H I1 + H b , and the interaction V be- 
tween them. In the last step of Eq. (3), the Liouville repre- 
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sentation in which the total Liouville operator L corresponds 
to the absorber molecule, the bath, and the interaction 

L=L a + L b + L l =L 0 +L\ (4) 

has been introduced. In terms of the Liouville representation, 
the spectral density can be written as 

F(w) = - Re Tr f V'(av*‘‘ X 'p/0^- (5) 

IT JO 


In the above expression, the superscript ( H ) of the Heisen- 
berg representation and the zero-time argument of /4 H, (0) 
have been dropped. It is well known that since the isotropic 
interaction does not depend on the rotational quantum num- 
bers, the Liouville operator associated with it will cancel out 
when it acts on the line space. Therefore, the isotropic inter- 
action need not be included in L, . 

We assume that 


-iLl s . e ~iL 0 l e -iL,l 


( 6 ) 


The limitation of this approximation has been discussed in 
previous papers. 7,11 Then F(w) can be expressed as 

1 

F (a))~- Re Tr {/ i a ■ exp[i(a>-L 0 )t]e~‘ L '‘p/t a }dt 

TT Jo 

1 1 
=-Im — t 
V 2 TTl J — 


xTr fi a 


dco' 

1 


—LA <d' -L 


P (a V, 


where 


cu'-L, 


= Td 


cu'-L, 




(7) 


( 8 ) 


and the summation indices of the trace a and b indicate the 
absorber and bath molecular variables, respectively. In the 
above expression for F(oj), at least formally, the bath vari- 
ables have been isolated. The resolvent operator l/(cu — Z-j) 
is defined only in the representation constructed from its 
eigenvectors. Therefore, we have to find the basis in which 
L i is diagonal. In general, the anisotropic interaction consists 
of several components and each of them has a specific 
spherical symmetry and r dependence. This implies that the 
diagonalization would have to be carried out for every r, but 
this is not feasible in practice. Therefore, in previous 
work 1-13 an assumption has been made that only one com- 
ponent of anisotropic interaction is dominant and all others 
can be ignored in the calculations. 

In order to overcome this fundamental drawback, a new 
method is introduced in the present work in which several 
components of the anisotropic interaction can be taken into 
account. We assume that the Liouville operator L t corre- 
sponding to the anisotropic interaction can be expressed as 
the sum of several terms with different spherical symmetries 
and r dependencies 

L i =L\ 0) + L\' ) +--- = GR l (r) + FR 2 (r) + --- (9) 


For clarity, only two terms are included explicitly in L) in 
the following expressions, although the inclusion of more 
terms is straightforward. In general, L\ 0) and L\ l) represent 
the leading and the next leading terms of the anisotropic 
interaction, respectively. Consequently, we have 


GR t (r)>FR 2 (r). 


( 10 ) 


We denote by and G ae ^G a -G 0 ) the 

eigenvectors and eigenvalues, respectively, of the Liouville 
operator G in line space notation constructed from | a) and 
G a , the eigenvectors and eigenvalues of the related matrix 
operator G, i.e., 


G\ a f) = G al3 \ a {3'). 


(ID 


Similarly, for the Liouville operator F, we have 


^lys^/vlys 7 ), 


(12) 


where | , y5 t )( s |'y)(4) F yS (^F y -F s ) are the eigenvec- 
tors and eigenvalues, respectively, of the Liouville operator F 
constructed from | y) and F y , the eigenvectors and eigenval- 
ues of the related matrix operator F. 

We make the approximation 


g~‘Llt~g it j 'e - '^ 




(13) 


which is valid if the contributions from the commutator 
} ] are small. It is worth mentioning that the cutoff 
on the overlap factor introduced later is consistent with this 
approximation. Introducing both Laplace transformations of 
e i and e , we have 


( Jrrij / J d<o" exp[- i(co' + u)")t] 


1 1 

X— nsrX- 




( 14 ) 


Then, the spectral density F( co) can be expressed as 
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xf Tr 
Jo 


Ma ■ exp[ i ( 01 - a ' - 01 " - L 0 ) t ] 




1 / l \ 2 r® r- 

= Im — dcj'l dm" 

tr \ 2m J - * J- co 


XTr 


Ma - 




‘ J '-1 

7r l 2-n-f J 


= Im 


1 /* 00 /*oo 


dee" 


X Tr ii a 


(J)— 0)' — (1)” — L n 




(15) 


Thus, we obtain an expression of F(oi) in which the bath 
variables have also been isolated analogous to Eq. (7). 


B. The matrix elements of the resolvent Liouville 
operator ( 1 /( a> ' - L H x 1 /( m" - L V 1 >)) 6 

The resolvent operator (l/(co’ -L^)x l/(a/'-L ( i 1) )) i> 
in Eq. (15) is an operator acting on the line space of the 
absorber molecule a. Therefore, in general, it can be ex- 
pressed as 


1 


1 


(o'-L\ 0)X w"-L\ [) l b 

Hi i Vi ' 1 1 ' b 


(16) 


By explicitly writing the summation over the bath molecular 
variables, its components can be written as 

= 2 < (< l7 l )(, 2 , 2 ) l^7 — /,(0> X pib) 


We note that the summation over the bath variables in line 
space in Eq. (17) includes only the states |i 2 i'|) because of 
the assumption that the bath molecule is not involved in the 
transition. 

In the present work, the translational motion is not 
treated quantum mechanically; rather one uses the quasistatic 
approximation. The classical ensemble average over the 
translational morion indicated by a subscript r is represented 
by an integration over r, the separation between two inter- 
acting molecules, with the statistical weight 
n b exp[- VanisJryfcTjexpf- V- lso (r)/kT] for a pair of mol- 
ecules, where n b is the number density of the bath molecules. 
Since the anisotropic potential consists of two terms, each 
contributes to the total statistical weight which can be ap- 
proximately expressed as a product form 

expC-V^iW] 

— exp[ - V^(r)/*r]exp[ - V^Jr)/*?’]. (i 8) 

In line space, the weight factors corresponding to each com- 
ponent of the anisotropic interaction can be expressed as op- 
erators constructed from their eigenvalues and eigenvectors 

ex P[“ v lSL( r ) /fcr ] = 2 

a$ 

(19) 
and 

exp[ - Vi^rykT] = 2 | y^)exp[ - F^rJ/kT^ y8 <\ . 

yS 

(20) 

Therefore, the ensemble-averaged component of the resol- 
vent operator (l/(«' -L ( , 0) )x \!{oi" -L\ u )) b r can be ap- 
proximately written as 


= 4^71*2 (O'lyiX'VDlf dr r 2 exp[- V lso (r)/kT] 
‘ii2 Jo 

—r~rm ex p[ “ KL(ry kT ] 

CO Lj 

x "^ZTTo ex PC — v Uiso( j P (&) l ( * [7 1 + ) O 2 / 2 ) ) 

= 47rn i ,2 2 E {{i\j\)(hii)\ a ^)f dr 7-2 


*2/2 a & yS 
X exp{ -[V iso (r)+ V , anjso ( ay;r)]/kT} 
1 1 


xKhj'SJUzjl))- 


(17) 
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lim =P — iir8(<o-L), 

c—Q id — L + ie (o — L 


(23) 


In the above expression, the factors (a)3 t j')'5 t )(=(aj7X^i)3)) 
describe the overlap between the eigenvectors of the Liou- 
ville operator G and the eigenvectors of the Liouville opera- 
tor F. If G and F commutate, they have common eigenvec- 
tors and {a(£\ y<5 f ) = 8 ay S^ & . 

To obtain an expression for Im{iiy|j{l/(o*' — Z.®) where P indicates the Cauchy principal value of the integral 

X \/(co" — we use the well-known formal and <5 is a Dirac function with a Liouville operator as its 


identity 


argument, and find 


Im< '’ iy " ‘ 1 \ w ' - If X - L‘ 1 V b ‘ f) 

= — 4ir 2 n b 222 <(‘Ji)(‘2‘l)! a! 0 t )[ dr r 2 exp{ - [ V lso ( r ) + ay,r)]lkT\ 

i, i, nr R * 0 


X \P 


12/2 “0 r s 
1 


'-G a 0 R { {r) ^ F ^ G °( >Rll ' r ^ P oj"-F YS R 2 ( r ) 


(a/8 t |y5 t )<y5 t |p (4) l( | 'u! t )0'2y2)) 


=-S S S <(»‘i./i)0V2)M t > 

1*2/2 a P y g 


o> l -G a f S R l [R^{< t >"IF yS )] \F yS \ 


o-r w 
R 2 | F 


Xexp 

XP 


-v., 


ay;R 


-i 


co 


2 \ F 


yS 


kT + 


\G a0 \ 


CO 


r a0. 


exp 


— V, 


yS 


-ii w 

«r;^i br- 
G a 0 


kT 


<D"-F rS R 2 [R; l (u>'/G af} )] 




(24) 


where R 2 1 (<d"/F 7S ) is the positive solution of the equation co"-F yS R 2 (r) = 0 and R { '(w'/G^) is the positive solution of 
the equation id' -G a pR { {r) = 0, respectively. In the above expression, u(r) and v(r) are defined by 


u(r) = 4v 2 n b Q(r) exp[ - V iso (r)/kT] 


(25) 


and 


u(r) = 4ir 2 /i 6 9(r) -j exp[- V lso (r)/kT], (26) 

where Z? '(r) = dR(r)/dr and O(r) is the unit step function defined as follows: 

©(r)= 1, if r> 0, 

6(r) = 0, if r< 0. 

C. The individual line shape functions 

To simplify the notation, we omit the subscript 1 of i, and j j indicating the absorber molecule. With Eqs. (15) and (24), 
the spectral density can be rewritten as 


F(w)= ~i 2 2 OlAli)^) f_ a 

Li W l,i 


= ^2 2 (/l/» a |*)-2 2 2 J ^ dcD'J ^dat" 


1/ 1 / 


<2h a & yS 


<D— (D —(D —ID: 
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X P 


(«*/f*)] \F yS \ 


R 2\ F 


0 ) 


y& 


|C 


a$\ 




(X) 


’« 0 


expj ^anis 




01 


r *0 


exp 
/ kT 


— V . 

9 antsol 


ay;R 2 1 


r« 


kT 


XP 


o> n -F yS R 2 [R; l (w'/G a0 )] 


X(^| r3 t )(r5 tj p ( 4 ) j(r/ t )(ji/ -t )){i . lp (a)^ |/) _ 


( 27 ) 


By carrying out the integrations in Eq. (27) (see the Appen- 
dix), we obtain 

F(<o)= ^2 2 0lj“al0-2 2 2 ((»y t )(‘2*2)l«j9 t > 

O' <’/' ‘iJi <*P rs 

X / 0ys(<*>- <*>//)( I yS 1 ) 
x<y5 t |p ( * ) |(iV' t )(/2/5))<i'IP w Ai'U'> 

= 01/*al«)*2 2 2 (('7 t )(‘2‘2)l Q! ) 3t > 

0 i 2 <*P 

XI a0yS (o)- a)ij){ap\ y5 t ><7|p w p (a) /t a |^) 

= E S 2 2 «“l''2)< 1 Vm|y)0'*2ly3»* 

ij i 2 yS 

(28) 

In the last expression for F(m), the summation indices i and 
j exclude their magnetic quanmm numbers since the latter 
are included in the summation over the magnetic quantum 
numbers indicated by {m}. If we introduce the individual 
line shape function <a t j) defined by 

Xij(oj-aj u ) = (w- Wy) 2 2 222 ((«l»2)<i|/4m|y) 

l" 1 ! i 2 af$ y& 

X(jh\P))*laliy6(u- Uij)(<X&\yti < ) 

x(y|p ( V a VJ^/p,KI 2 ’ (29) 

where the reduced matrix elements p.y is defined by 

{i\p-m\j)~C{j\i-,m j mm i )(i ij l{2i+\) u ' 1 , (30) 

then the spectral density takes a more familiar form 

F(fc>) = ^ 2 PilMol 2 * y(«-»y). OD 


We can also introduce the symmetric individual line 
shape functions which correspond to the replacement 
ail a0yS (<i>- o>ij) by their symmetric partners I a07 s(u— <u fy ) 
(see the Appendix) in the expression for * ly (a>— a> /; ), 

*/,(<*>- a», v ) = exp[fi( w- <u y )/ 2 *r]xy(a>- a>y) 

= («>-<*>, y) 2 2 2 2 2 ((®l**2>(»|y* m |y) 

( m ( i 2 ap yS 

x 0‘«2!i8»*/^ r «(a»- tVyXa^l r<5 f ) 

X(7|p< 4) p (o V m |<5)/p,|/iy| 2 . (32) 

To satisfy the detailed balance requirement, 14 the spectral 
density must take the following form: 


F{(o)-— 2 exp[-ft(o>— <t)ij)/2kT] 






l 


(a >+(i>ij) 


jXji(<t> + Wij) 


Pi\Pij\ 


with 


X'j(°>) = Xji(-u)- 


(33) 


(34) 


Consequendy, in terms of ^y(cu), the total absorption coef- 
ficient a(<v) is given by 


a{o)) = n a 2 S>j 

*ij> 0 


< osh(h(o/2kT ) 1 

1 

uiijSh(ha>;j/2kT) it 

(co-Uij ) 2 


X Xij(<o- o»y) + (<a + Wy) j */,(<*>+ «y) 
where 5 i; - is the usual line strength. 15 


(35) 


J 


D. The band-averaged line shape function 

For practical calculation, we introduce the intensity weighted band-averaged line shape function defined by 


t • 

* m -2 Pi\Pij\ 2 Xij(°>) j 2 PilPijl 1 

ii > ii 


= w 2 2 2 {a\p m \P)*i a0y6 M{a^\yS^){y\p w p {a) p m \S) /2 P/I/4, yl 2 

a 0 yS • U 
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= 4'JT 2 n 6 0> 2 2 X ( a l/ A ml^>*©( r c) 

a0 yS 


1 

G a pR\{r c ) + F yS R' 2 {r c ) 


r 


2 

c 


X «p{ - V ao {r c VkT-[{G a +Gp)R i {r c ) + (F y + F Mr c )]/2kT}(afi t l yS0< y|p ( V° V J <5> 



(36) 


r 


where r c is determined by 

<o- G afi R,(r c ) + F yS R 2 (r c ), 


(37) 


or approximately by 


r c =R7' 




FysR 2 




(38) 


In general, for a specified pair of a and /J, the value 
of r c and the weighting factor exp{- 
IkT- [(G a + Gp)R { (r c ) + (F y + F s )R r (r c )]/2kT} depends 
on y and S. However, in practice, due to the overlap factor 
(a/3 t | y<5 7 ), the summations over y and 8 can be limited to 
those pairs which significantly overlap with a and /? only. It 
turns out that the most overlapping pair of y and 8 is that 
whose eigenvalues F y and F s are the closest to 
F „^(a\F\a) and F p^(j3\F\(3), respectively. In fact, we 
can express F a and Fp as 


X{o>)=47T 2 n b ai 2 'Z HaU/ 3>*©(7 c ) 

ct$ 


\C^\M + F^U\ f '”*- VUf ’ VkT 

- [( G a +G p )R x Cr c ) + {F a + F 0 )R 2 (r c ) ]/2 kT} 


X<a|p<VVj/3> 



(41) 


More suitably, we can also introduce the positive and 
negative frequency resonance-average line shape functions 
£ + (<u) and *_(cu) as we did previously. 7 We define the inten- 
sity weighted positive frequency resonance-average line 
shape function by 

*+M= X Xij{u)Pi\ni)\ z j X PiK'Vm |y)l 2 . (42) 

a>ij>0 / ij 

Ei>Ej 


<«l y)Fj{y\a) 
y 


(39a) 


and the intensity weighted negative frequency resonance- 
average line shape function by 


and 


*>-X (mpsim. (39b) 

s 

respectively. Since they are the average values of F y and F s 
over y and 8 with the overlapping weights |(a|y)| 2 and 
|(^/3)| 2 , they are more likely close to the most heavily 
weighted values. The more overlapping of {a0\y S f ), the 
closer is the average values and the most heavily weighted 
values. Then as an approximation we can replace F y and F s 
by F a and Fp , and replace (a/3 f |y<5 f ) by ^ayfys ■ These re- 
placements can also be made for the case of 
GR [ (r)'i>FR 2 (r) since the usual perturbation method is ap- 
plicable. In these cases, the corresponding value of r c de- 
noted by r c , 


F C = R 


- 1 
i 




F a 0 R i 


R 



(40) 


and the weighting factors are independent of y and 8. Then 
the expression for x(oj) can be simplified by omitting the 
summations over y and 8, i.e.. 


X-M- X XjiMPj\p-ji\ Z / X PylOlAml')! 2 - (43) 

/ ij 

E,>E, 

Their explicit expressions are easily obtained by comparing 
with the expression for x(io), but are not presented here. 
With the two average line shape functions, the absorption 
coefficient can be written as 


v ^ u)sh(h<i)/2kT) 1 1 

S ‘ J UijSh^ijUkT) n [(*>-*>,/ 


X^ + ((u-w, v ) + 


1 

(a > + cj u ) 


T x-(v+Uij) 


(44) 


If one does not distinguish between these two different func- 
tions, then in terms of the band-average line shape function 
x(<d), the absorption coefficient a(cu) is given by 


a>sh{hu)/2kT ) 1 

a(w) = /j a 2j T77 TTfyF ~ 

\ >o a>ijsh{fi(»,j/2kT) tt 


l 

(w - WijY 


Xx{o>-u> u ) + 


I 

(cO+OJij) 


7 i(0)+U >,;) 


(45) 
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III. SOME EXAMPLES 

As a first example, we consider a system in which the 
leading anisotropic interaction is the dipole^dipole interac- 
tion and the next leading term is the dipole-quadrupole in- 
teraction. This implies that 


r o= 



1/3 


(48) 


It is obvious that r 0 is the lowest-order approximation for r c . 
Then Eq. (A 12) can be expressed as 


R i( r ) = p an d R 2 (r)=-p; . 


(46) 


With the assumption G( l/r 3 )>f ( 1/r 4 ), it is easy to find the 
approximate values for r c in Eq. (38) 




1+ 3 \G 


a0, 


4/3 


FySU 


1/3 


(47) 


where 


. 4ir 2 1 

I a f3ys(u) = —nb 


2 

1+ 3 


FySU 


1 

G a 0 

X0(r f )|G^|exp{-V iso (r f )/^7 


4/3 


1/3 


- UG a + G 0 )/rl + (F y + F s )/r* c ]/2kT}. (49) 

Accordingly, we find the symmetric band-averaged line 
shape function 


4 7 j 

XM = —n b Z {cc\fi m \PrQ(r c )\G ap \2 tM~V^r c )lkT-[{G a + G p)!r] + (F y + F s )lr*]l2kT}{ccp\y8') 

*P yS 

{ r |p( i,) p (£,) / i m |5) + - |— ) <o' n (y\F yS p (b) p ia) jj. m \8) j X P/l/^oi 2 


4 77 

(ct\^ m \l3)*Q{? c )\G aP \txp{-VJ < ? c )lkT-[{G a +Gp)l?] + {F a +F 0 )l^ c \l2kT} 

a0 


/ X P.lM.yl 2 . 
/ ij 


(50) 


2 / 1 \ 4/3 

(a\p (b) p {a) p. m \/3)+- cu 1/3 (a|[F,p <i) p (fl Vm]l^> 

where F is the matrix operator and 

[F,p (6) p (a Vm] = Fp (4) p (a - p (i,) p ( ° V m F. 

In the last expression of Eq. (50), a simplified form is given by using average values. The lowest-order approximation to £(o>) 
takes a form corresponding to the one component case of the anisotropic interaction (the leading dipole-dipole term) consid- 
ered previously 7 

4-7 r 2 


4 77 

X (0, («)— (a\pJ/3)*Q{r o )\Gjexp{-V iso (r o )/kT-[(G a +G 0 )/rl]/2kT} 
«P 

X(a\p ib) p {a) p. m \/3) I Z P/lM/yl 2 - 
• U 


(51) 


As a second example, we consider a case in which the two components of the anisotropic interaction have the same r 
dependence, i.e., R t (r) = R 2 (r)=R(r). Then Eq. (A12) can be expressed as 


/ «Pr«( w ) = 4 ' n ' 2,, * e ( r c) 
1 


( G af3 +F yS )R'(r c ) 


r 2 exp[ - V,J r c )/kT—(G a +Gp+F r + F $)R(r c )/2kT] 


where 


r c = R' 


= ^b O' 9(r c ) 


CO 


R {r c ) 


R'(r c ) 


r] exp[ - V iso (r c )/kT-( G a + G„+ F y + F s )R{r c )!2kT] 


(52) 


G a p +F y5l 


(53) 


Then, 
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^(w) = 4'jr 2 n i jo>iX ©(r c ) 


R(r c ) 


I R’(r e ) 

x (a/S t |r-5 f )(rlp <A) p (a Vml^} /2 p.l^.jl 


'1 exp{ - V,J r c )lkT-(G a + G 0 +F y + F s )R{r c )l2kT} 


= 47r 2 n 2 ,|w|2 (a\fj. m \0)*Q{? c ) 




R{f e ) 


R'ifc) 


fe*P{ ~ Vise ( r c )lkT-{G a + G p+ F a + F p)R(r c )/2kT} 


X(a|p ( *V‘ J Vml^) / 2Pil/*/y| 2 . 


(54) 


r 


where 



cj 

G*0+F af}i 


(55) 


These expressions are consistent with those obtained directly 
from one component formalism. 

Finally, we consider a system consisting of a linear ab- 
sorber molecule and a perturber atom. The isotropic potential 
is assumed to be a Lennard-Jones model 


V iso (r) = 4e 




(56) 


where e and cr are the usual parameters. In general, the an- 
isotropic potential can be expressed as 


V , ani so(T,^) = 2 Plicos ff)V,(r), (57) 

i 

where 0 is the angle between the absorber-perturber axis and 
the linear molecular axis, and P t { cos 6) is the associated 
Legendre polynomial. For a highly anisotropic system, sev- 
eral components /’((cos ff) have to be included in the expan- 
sion of Eq. (57). As an example, we consider the C0 2 -Ar 
system where only even / occur because of the symmetry of 
C0 2 . Therefore, the anisotropic interaction takes the form 

^aniso(^> 8) — P 2 (cos 6) V 2 {r) + P 4 (cos 0)V 4 (r) 

+ P 6 (cos 0)V' 6 (r)+--- . (58) 

For convenience, we model V,0) by 


V|(r)- 4e 



12 

-A 



(59) 


where B t and A, are the parameters. Then the potential en- 
ergy surface of the system is simply determined by e, cr, B t , 
and A;. 

Taking the absorber-perturber axis along the z direction 
of a space-fixed frame, the orientation of C0 2 with respect to 
the space-fixed frame is represented by 9 and <p, and the 
anisotropic potential is independent of <p. Writing 
P ((cos 0)=[(4 tt)/(2/+ 1)J Y !0 ( ff,<p), it is easy to diagonize 
the components P ,( cos 9) in the rotational basis |Im) since 
all the matrices of the components /’((cos ff) take exactly the 
same block form, and each of these subblocks is labeled by 
the magnetic quantum number m. Then, in the subspace la- 
beled by m, one can use the notation \am), | ym), ]fm),... 
and G a , F y , Z/ fl ... to represent the eigenvectors and the 
eigenvalues of the components P 2 ( cos ff), P 4 (cos ff), 
P 6 ( cos ff), etc., 

P 2 (cos ff)\ocm) = G a \am), 

P 4 (cos 9)\ ym) — F y \ ym), (60) 

P (,{ cos 9)\Cm) = H ( \£m). 

For simplifying the notation, the label m will not be written 
explicitly in further expressions. Including the three compo- 
nents given above, one is able to derive an expression for 
1 a0yS(v {oj- similar to Eq. (A12), viz. 


hpySfvi w ) = 4 •7r 2 n*0( r c ) 


1 

G a0 V’(r c ) + F y5 V:(r c ) + H iv V^r c ) 


r 


2 

c 


Xexp{-V iso (r c )/kr-[(G a + G^)K 2 (r c ) + (F r +F J )V 4 (r c ) + (// i +//,)K 6 (r c )]/2kr}, 


where r c is determined by 


(61) 


u = G a0 V 2 (r c ) + F yS V 4 (r c )+H iv V 6 (r c ). 


(62) 


Then, one can write the symmetric band-averaged line shape function as 
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FIG. 1. The total potential energy surface (in units of Kelvin) for C0 2 -Ar [Eq. (70)] given by Boissoies et aL (Ref. 13). The distances are in atomic units and 
the angles are in degrees. 


X(<*>) = 4 tt 2 n b a> 2 ^ (a\p . m \ £)*2 2 0(r c ) 

a/3 yS {tj 


G 0 ^2(r c ) + F r5 v;(r c ) + // f ,v;(r c ) 


X exp{ - V lso (r c )/kT- [(G a + C 0 ) V 2 (r c ) + (F r + F s ) V 4 (r c ) + 4 - HJ V 6 (r c )]/2kT} 

X(a/3 1 \yS t )(yS t \CT] t )(C\p {b) p {a) fi m \v) / 2 Pi\Pij\ 2 - 


(63) 


In the above expression, (a$\y8 x ) and {y<5'|£V) are the overlap factors. As we mentioned previously, for a specified pair of 
a and /?, the summations over the pair of y and S and over the pair of £ and rj can be limited to those which are most heavily 
overlapped with the pair of a and (3 only. In the case of V 2 (r)^>V 6 (r) and K 4 (r)> V 6 (r), one can simply set £=a and 7 =/3 
and replace r c which depends on a, /3, 7 , 8, £, and 77 by their approximated values r c determined by 

a>= G a pV 2 (r c ) + F r jV 4 (f c ) + H a pV 6 (r c ), (64) 

where 

H a0 =(a\P 6 (cos 0)\a)-(/3\P 6 {cos d)\/3). (65) 

This implies that r c is independent of £ and rj. Consequently, the symmetric band-averaged line shape function can be 
approximated by 
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£(a>) = 47r 2 n b u) 2 ^ <a|/rJ)3)*(a|p ( *V (a) /i m |/3)2 (a)8 t |75 t ><75 t |a/3 t )e(r c ) 

a/3 Y S 


G afi VZr e ) + F yS V 4 (r e ) + H a0 V 6 {r e ) | ‘ 

Xexp{- VJLr c )/kT-l(G a +GJVi(r c ) + (F r +FJV 4 (r c ) + {H a +H fi )V 6 {r c )y2kTl / 2 P,Im, 7 | 2 - 

/ I# 


( 66 ) 


Furthermore, if V , 4 (r)> K 6 (r), one can set y=a: and £=/3 further and replace r c by the approximate values r c 

determined by the following equation: 


<o = G a0 V 2 (r c ) + F a/3 V 4 (r c ) + H a0 V 6 { F c ) , (67) 

where 

^ap=(«l^ 4 (cos 0)|or>-(0|F 4 (cos 6)\f3). (68) 

Then, r c is independent of y, 8, £, and rj. Consequently, the symmetric band-averaged line shape function can be simplified to 



FTG. 2. The anisotropic pan of the potential energy surface (in units of Kelvin) for C0 2 -Ar [Eq. (70)] given by Boissoles et al. (Ref. 13). The distances are 
in atomic units and the angles are in degrees. 
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FIG. 3. The band-averaged line shape function (in units of cm' 1 atm' 1 ) as 
a function of frequency to (in units of cm" 1 ) for lines in the Ar-brdadened vy 
band of C0 2 calculated for 7=296 K using the potential surface shown in 
Fig. 1. 



FIG. 4. A comparison between the calculated absorption coefficient (in units 
of cm 2 molecule' 1 atm' 1 ) for 7=296 K in the 2400-2580 cm' 1 spectral 
region of C0 2 -Ar with experimental data. The solid line is the theoretical 
values calculated by Boissoles et al. (Ref. 13) using the potential energy 
surface given by Eq. (70); the dashed curve is the results obtained assuming 
Lorentzian line shapes; and the experimental results (Ref. 13) are denoted 
by +. 



FIG. 5. The function V 0 (r) (solid line) vs r obtained by fining a Lennard- 
Jones model to the data of Parker and co-workers (Ref. 17) (denoted by +). 
All quantities are in atomic units. 



FIG. 6. The same as Fig. 5 for the function V 2 (r). 
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*M=4 it 2 n b u > 2 2 

a0 




G a0 V’ 2 (F c ) + F a pV' 4 ( r c ) + H a0 V' 6 ( F c )\ c 

Xtxp{-V iso (F c )/kT-[(G a +G 0 )V 2 (F c ) + (F a +F 0 )V 4 (F c ) + (H a + H 0 )V (> (F c )]/2kT} /]£ p,|/t tf | 2 . (69) 

' ij 


IV. APPLICATION TO THE v 2 BAND OF C0 2 
BROADENED BY Ar 


We consider the v 3 band of C0 2 broadened by Ar in 
which both the experimental measurement of the absorption 
coefficient and some knowledge of the interaction potential 
between C0 2 and Ar are available at the present time. The 
wing of the Ar-broadened v 2 band of C0 2 for 7=296 K has 
been studied by Boissoles et a/. 13 recently using an effective 
“one-component” potential determined empirically from the 
spectral moments of the infrared vibration-rotation bands. 
The form of the potential used in their calculations is 


V(r, 8) = 4e 



6 


+ 4e\ 


cr\ 12 / <r' v6 

B 2\-j -^7 


Pli cos 8), (70) 


where the parameters are e/k= 153 K, <r=3.9t A, 
B 2 =l.09, and A 2 = 0.266. 

We note that the formula used by Boissoles eta/. 13 to 
calculate the band-averaged line shape function is similar to 
our formula obtained previously corresponding to the one- 
component anisotropic potential case. However, they have 


introduced an energetically corrected function in their for- 
malism. As we discuss elsewhere, 16 this correction corre- 
sponds to rewriting the spectral density in a more symmetric 
from 


F(to) = -ReTr f e ,m [n a 4p-e~‘ Lt Vp nMt- (71) 
ir Jo 

Similarly, we follow the same energetically corrected proce- 
dure in the present calculations, but we do not reproduce the 
expressions here since they are easily obtained from the ex- 
pressions given above. For the frequency of interest here 
2400-2580 cm -1 , we consider the contributions from the v 3 
band of C0 2 only and exclude any contributions from other 
bands. From the HITRAN data base, 15 there are 109 line listed 
in this band and all of them are included in the calculations. 
In Figs. 1 and 2, we present the total potential energy surface 
and the anisotropic potential energy surface of the C0 2 -Ar 
system, respectively, given by Eq. (70). In Fig. 3, the theo- 
retically calculated band-averaged line shape function for 
T- 296 K based on the potential used by Boissoles et al. is 
shown. In Fig. 4, the comparison between the calculated ab- 
sorption coefficient and the experimental data for T= 296 K 




FIG. 7. The same as Fig. 5 for the function V 4 (r). FIG. 8. The same as Fig. 5 for the function V' 6 (r). 
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FIG. 9. The total potential energy surface (in units of Kelvin) for COj-Ar obtained by fitting a Lennard-Jones model to the data of Parker and co-workers 
(Ref. 17). The distances are in atomic units and the angles are in degrees. 


is presented. The results obtained from a Lorentzian line 
shape are also given in the figure. 

For the potential energy surface of the CCL-Ar system, 
there are some theoretical calculations done by Parker 
et at . 17 In their calculations, the short range potential consists 
of two terms — V HF (r,d) to take into account the Coulomb, 
kinetic, and exchange contributions to the interaction energy 
and V COR (r,0) to take into account the correlation contribu- 
tion. The values of V HF (r, d) and V COR (r,0) obtained are 
tabled as the functions of distance and angle (see Tables III 
and IV of Ref. 17). By fitting these results with an expansion 
of the Legendre polynomials at each r, 

Vhf ('•.*)= 2 V? F (r)P,( cos e ), 

1 = 0 . 2 . 4 .... 

(72) 

V C qr ( r.0)= 2 Vf OR (r)P,(cos 9), 

/= 0 , 2 . 4 .... 

the values of V ( HF (r) and vf 0 R (r) at this r are obtained and 
listed in Tables VII and VIII respectively, of Ref. 17. They 
have further modified the values of vf 0 R (r) obtained to take 
into account the proper van der Waals behavior at large dis- 
tances. However, since the far-wing line shapes are mainly 
determined by the short range behavior of the interaction 


potential, we will not follow their next step and only use the 
values given in these tables in our calculation. 

For simplifying our calculation, analytical expressions 
for Vj^ir) and vf 0R {r) instead of the tabulated values are 
necessary. In addition, for later convenience, we prefer to 
have them take the form of the generalized Lennard-Jones 
model mentioned previously. It turns out that one is able to 
fit some coefficients pretty well and to fit others only reason- 
ably well by adjusting the parameters of the Lennard-Jones 
model. We note that there are some uncertainties in the pa- 
rameters obtained due to the difficulty of the fitting. In the 
present calculation, we consider the terms up to 1 = 6 only 
since the magnitudes of the higher terms are much smaller 
than the lower terms. In Figs. 5-8, we plot the results of 
fitting for V 0 (r), V 2 (r), V 4 (r), and V 6 (r). The potential pa- 
rameters obtained are elk = 62.0 K, cr=4.44 A, B 2 — 1.44, 
A 2 = 0.73, S 4 = 0.70, ^4 = 0. 13, B 6 = 0.19, andA 6 = 0.03. 
In Figs. 9 and 10, the total potential energy surface and the 
anisotropic potential energy surface corresponding to the po- 
tential parameters given above are presented. Based on this 
potential, we first calculate the energetically corrected band- 
averaged line shape function for 7=296 K from Eq. (69) in 
which the summations over the pair of y and 8 and the pair 
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FIG. 10. The anisotropic part of the energy surface (in units of Kelvin) for C0 2 -Ar obtained by fitting a Lennard-Jones model to the data of Parker and 
co-workers (Ref. 17). The distances are in atomic units and the angles are in degrees. 


of f and 7] have been averaged out and only the summation 
over the pair of a and /3 is performed. Then we repeat the 
same calculation with Eq. (66), in which only the summation 
over the pair of f and rj has been averaged out and the 
summation over the pair of y and S is performed correspond- 
ing to those pairs with an overlapping restriction, i.e., a cut- 
off. This cutoff is necessary to be introduced in calculations. 
It represents a measure of how wide the averaged angular 
distribution of the C0 2 wave functions with the same mag- 
netic quantum number m spreads. However, in the present 
study, we do not discuss it further and simply choose 
|(a|y)| 2 X).01 and |<<5)/3>| 2 ^0.01 . Since the eigenvectors of 
/ > 2 (cos 6) more or less likely overlap with the eigenvectors 
of P 4 (cos 6), for a specified |a), the summation over y with 
the restriction |{ar|y)| 2 3s0.01 includes at least 94% of the 
states \y). We present the line shape functions obtained in 
- Fig. 1 1 and the corresponding absorption coefficients in Fig. 
12. For comparison, the results obtained with the “one- 
component” model of Boissoles et al . 13 are also presented in 
Fig. 12. 

It is well known that the profile of the line shape func- 
tion is determined by the interaction potential surface. How- 


ever, the behavior of the far wing of the line shape function 
is mainly governed by the short range interaction. More spe- 
cifically, first of all, the anisotropic potential energy surface 
determines beyond what kind of positions labeled by r and d 
for the C0 2 -Ar system are able to contribute to the values of 
shape function at a specified tu. Second, the total potential 
energy surface determines how much will contribute to it 
from these locations. For example, in order to estimate the 
relative magnitude of the shape function around (o = 200 
cm" 1 in comparison with that around 50 cm -1 , one first has 
to locate the areas in the anisotropic energy surface in which 
the absolute values of anisotropic energy are at least equal to 
or larger than the corresponding u> values. Then, in a graph 
of the total energy surface, one has to find out whether there 
are significant probabilities to locate a pair of C0 2 and Ar at 
these areas. Therefore, it is not the total potential energy 
surface alone, but the relationship between the total and the 
anisotropic interaction energy surface to determine the pro- 
file of the far-wing line shape. These conclusions are easily 
understood since the first one comes from a basic require- 
ment of the energy conservation and the second one is re- 
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FIG. 11. A comparison between the line shape functions (in units of cm' 1 
atm' 1 ) for T = 296 K obtained with the multicomponent anisotropic model 
given by Eqs. (66) (solid line) and (69) (dotted line), respectively. 

lated to a statistically probability to find a pair of C0 2 and Ar 
at a specified location. 

V. DISCUSSION AND CONCLUSIONS 

As can been seen clearly from Figs. 3, 11, and 12 that 
there are substantial differences between the single- 
component and the multicomponent results for both the line 
shape and the absorption coefficient. The differences be- 
tween the results obtained from the two different multicom- 
ponent formulations are small. The present theory can repro- 
duce accurately the experimental data for the absorption of 
the high-frequency wing of the v 3 band of C0 2 broadened by 
Ar. Because this spectral region lies beyond the bandhead, 
most of the absorption arises from the far wings of lines. 
Therefore, in the present calculations, we have not included 
either the near- wing corrections" resulting from the non- 
commutation of the Liouville operator describing the unper- 
turbed system and that describing the interaction [cf. Eq. 
(6)], or corrections arising from molecular motion. 12 Of 
course, for the accurate calculation of the near wing of the 
lines or absorption resulting from this region (e.g., near the 
bandhead or closer to the band center), these refinements 
should be included. 

In the present work, we have introduced a cutoff for the 
overlap factors and arbitrarily chose the values 0.01 for all 
the subblocks labeled by m. In fact, the cutoff depends on m 
since the cutoff is related to the averaged angular distribution 
functions of the C0 2 wave functions with the same magnetic 
quantum numbers. On the other hand, we have fit the tabular 
data for the C0 2 -Ar interaction to a Lennard-Jones model to 
facilitate the calculations. However, this model lacks the 



FIG. 12. A comparison between the calculated absorption coefficient (in 
units of cm 2 molecule' 1 atm' 1 ) for T= 296 K in the 2400-2580 cm' 1 
spectral region of C0 2 -Ar with experimental data. The solid line is the 
results obtained using Eq. (66). the dotted line is the results obtained using 
Eq. (69), and the dashed line is the results of Boissoles el al. (Ref. 13). 

flexibility to present the r dependence well especially for the 
low-order (e.g., I = 0 and 2) components. Thus, there is some 
uncertainty introduced by this procedure as well as in the 
original data, and these uncertainues affect the magnitude of 
the absorption. In fact, the absorption is very sensitive to the 
details of the interaction potential. 

With the improvement in the theory detailed in the 
present paper, one is now able to relate the absorption in the 
far-wing region of pure rotational or vibration-rotational 
bands directly to the interaction potential. Although for sim- 
plicity in the present paper we considered a linear molecule 
perturbed by an atom, the basic theory is also applicable to 
more complicated problems, e.g., molecule-molecule sys- 
tems. We are in the process of making such calculations, and 
after the results are compared with experimental data, we 
will be in a better position to assess the overall accuracy of 
the theory. Nevertheless, we feel that the present theory pro- 
vides a consistent framework for the calculation of the far 
wings of spectral lines within the binary collision approxi- 
mation. Together with improved impact models valid near 
the centers of lines, we can now calculate accurately the line 
shape for both small and large frequency displacements from 
the line centers. However, there is still considerable uncer- 
tainty for intermediate displacements that must be resolved 
before one can attain the goal of computing accurate line 
shapes over their complete extent. 
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appendix 

In order to carry out the integrations denoted by 
lafiy ^ Eq. (27), 
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we have to keep in mind that both the integral variables co' and <o" have small imaginary parts icr'(a'>0) and icr”(cr"> 0) 
which have not been expressed explicitly in the formalism and that a' and <f will go to 0 + after the integrations are completed. 
Also we have to keep in mind that the frequency co has a small imaginary part ie (e>0) too and that finally e will go to 0 + after 
all of the calculations are completed. Now we begin to calculate the first term in Eq. (27), i.e., 
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(A2) 


We carry out the integration over co' first and construct a closed curve C which consists of three parts: (1) a straight line along 
axis icr' from -°° + icr' to G a gR t [R 2 \co"/F yS )]- S+ icr' and from G a gR i{R 2 x (co" ! F yS )] + 8+ ia' to » + icr'; (2) C s , a 
half-circle with the point G a pR x [R^ \co"/F yS )] + icr' as its center and with a radius 5 on the upper plane; (3) C R , a half-circle 
with the point icr' as its center and with an infinity radius on the upper plane. Then we have 
1 1 
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With Eq. (A3), Eq. (A2) can be written as 
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In order to carry out the second integration over co", we introduce a new integration variable y—y{co") defined by 
co 
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or alternatively given by the pair of equations 
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y = G a0 R 1 ( r c ) + F yS R 2 ( r c ) . 

We have 
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With Eqs. (A 5) and (A7), Eq. (A4) can be rewritten as 
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where r c is determined by 

a)-bJij=G a 0Ri(r c ) + F yS R 2 {r c ). (A9) 

We note that r c will be the same if the value of (cj-cj^) is changed by a sign and both of the indices a and /3 and the indices 
y and S are interchanged simultaneously. 

Similarly, we can calculate the second term in Eq. (27), i.e., 

(ai) r„ d “’L / kT . 


r u>"-F yS R z [R; \u>' I G afi )}-' V ' 

However, by comparing Eqs. (A2) and (A 10), one is able to find out that one integrand becomes another integrand simply by 
interchanging between one set of symbols to' co" a 0 y S F G R { R 2 and another set of symbols id" cd’ y8a/3GFR 2 R l , 
simultaneously. In addition, we note that the final result of the first part integration given in Eqs. (A8) and (A9) is independent 
of interchanging between one set of symbols a /3 y 8 F G R x R 2 and another set of symbols y 8 a j3 G F R 2 R x . Therefore, 
one concludes that the contributions from these two parts of integration are the same and is able to find 

^rs(“-Wi / ) = 4lrZ "i 0 ('-c) ^ ' p/7 ' vie p'T T r c ex P{“[ V i S o( r c) + G a F 1 (r ( .) + F r F 2 (r c )]/l:T}. (All) 

u a0 K lK r c> + f ySKzCc) 

For later convenience, we can introduce symmetric functions ~I a p y6 (u)— o>^) defined by 
hpydi u - a>ij) = exp[ft ( w - (Dij)/2kT]I a p yS { id - <a i} ) 

^^v{-y^{r c )lkT-[{G a + G e )R l (r c ) + {F y +F s )R 2 (r c )]l2kT}. (A12) 

It is worth mentioning that these functions have a symmetric property 

hpysi W - 0>ij) = 1 0aSy( ~ « + “,,) • ( A 1 3 ) 

Based on the knowledge about r c given previously, one is able to derive this property easily. 
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